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Abstract 

The diffractive dissociation of virtual photons, 7*p — ► Xp, has been studied 
with the HI and ZEUS detectors at HERA using various complementary tech- 
niques. Events have been selected by direct tagging of the outgoing proton 
or by requiring a large rapidity gap between the proton and the system X. 
The diffractive contribution has also been unfolded by decomposition of the 
inclusive hadronic final state invariant mass distribution. Here, detailed com- 
parisons are made between diffractive cross section measurements obtained 
from the different methods and the two experiments, showing them to be con- 
sistent within the large uncertainties associated with the treatment of proton 
dissociation processes. First steps are taken towards the combination of the 
HI and ZEUS results. 

1 Introduction 

In the single diffractive dissociation process in proton-proton scattering, 
pp — > Xp, at least one of the beam hadrons emerges intact from the col- 
lision, having lost only a small fraction of its energy and gained only a small 
transverse momentum. In the analogous process involving virtual photons, 
7*p — > Xp (figure 1) [2,3], an exchanged photon of virtuality Q 2 dissociates 
through its interaction with the proton at a squared four momentum transfer t 
to produce a hadronic system X with mass Mx- The fractional longitudinal 
momentum loss of the proton during the interaction is denoted xjp, while the 
fraction of this momentum carried by the struck quark is denoted j3. These 
variables are related to Bjorken x by x = (5 xjp. 

Diffractive interactions are often discussed in the framework of Regge 
phenomenology [4] in terms of the exchange of a 'pomeron' with vacuum 
quantum numbers. This interpretation in terms of a universal exchange is 
experimentally supported by the 'proton vertex factorisation', which holds 
to good approximation over much of the accessible kinematic range at low 
xjp, whereby the dependences on variables describing the soft interaction 

with the proton (xp, t) factorise from those related to the hard interaction with the virtual photon (/?, Q 2 ). 
Similar reactions, in which sub-leading Reggeon and pion trajectories are exchanged, have a negligible 
cross section at the smallest values. 

Significant progress has been made in understanding diffraction in terms of QCD by studying 
virtual photon dissociation in deep inelastic electron-proton (ep) scattering (diffractive DIS) at HERA. 
A recent review can be found at [5]. As well as being sensitive to novel features of partem dynamics 
in the high density, low x regime, diffractive DIS cross sections are used to extract diffractive partem 
density functions (DPDFs) [6-10], an essential ingredient in predicting many diffractive processes at the 
LHC and in estimating backgrounds to more exotic processes such as central exclusive Higgs production 
{pp -» pHp) [11]. 




Fig. 1 : Illustration of the kine- 
matic variables describing the 
virtual photon dissociation pro- 
cess, 7*p — > Xp, in ep colli- 



* Contributed to the Proceedings of the HERA-LHC Workshop [1]. 



Similarly to inclusive DIS, cross section measurements for the reaction ep — » eXp are convention- 



ally expressed in terms of the reduced diffractive cross section, a, 
cross section by 
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At moderate inelasticities y, cy*^ corresponds to the diffractive structure function to good ap- 

proximation. In this contribution, we tackle the technical issue of compatibility between different a® 
data sets through detailed comparisons between different measurements by the HI and ZEUS collabora- 
tions and take the first steps towards a combined HERA data set. 



2 Methods of selecting diffraction at HERA 

Experimentally, diffractive ep scattering is characterised by the presence of a leading proton in the final 
state retaining most of the initial state proton energy, and by a lack of hadronic activity in the forward 
(outgoing proton) direction, such that the system X is cleanly separated and Mx may be measured in the 
central detector components. These signatures have been widely exploited at HERA to select diffractive 
events by tagging the outgoing proton in the HI Forward Proton Spectrometer or the ZEUS Leading 
Proton Spectrometer (proton-tagging method [7, 12-15]) or by requiring the presence of a large gap in 
the rapidity distribution of hadronic final state particles in the forward region (LRG method [6,9, 16, 17]). 
In a third approach (Mx method [17-20]), the inclusive DIS sample is decomposed into diffractive and 
non-diffractive contributions based on their characteristic dependences on Mx- 

The kinematic coverages of the LRG and Mx methods are limited to xjp < 0.05 by the need 
to contain the system X in the central detector components. These two methods are equivalent for 
Mx — > 0, but differences are to be expected at larger Mx, where the LRG method measures the full cross 
section from all sources at a given (sp, (3, Q 2 ) point, whereas the Mx method involves the subtraction 
of a 'non-diffractive' component. LPS and FPS data extend to x-p ~ 0.1 and are therefore the most 
sensitive to non-leading contributions, including Reggeon and pion trajectory exchanges. Apart from 
the proton dissociation treatment in the HI case (see section 4.2), the cross sections measured by the 
proton-tagging and LRG methods are equivalent. 

The methods differ substantially in their dominant sources of systematic uncertainty. In the LRG 
and Mx methods, the largest uncertainties arise from the admixture of low mass leading baryon systems 
other than protons. These include proton excitations to low mass states as well as leading neutrons pro- 
duced via charge exchange reactions. All such contributions are collectively referred to here as 'proton 
dissociation', ep — > eXN, with the baryon state N having mass M^. Proton dissociation processes can- 
not always be distinguished by the LRG and Mx methods from events in which the proton is scattered 
elastically. Conversely, low-xp samples selected by the proton-tagging method have little or no proton 
dissociation background, but are subject to large uncertainties in the proton tagging efficiency, which is 
strongly dependent on the proton-beam optics. Proton spectrometers also allow a measurement of t, but 
the statistical precision is limited by their small acceptances. 

Comparing the results from the three different methods is a powerful test of the control over 
the systematics of the measurements. At low ip, the ratio of results obtained by the LRG and Mx 
methods to those from the proton-tagging method can also be used to quantify the proton dissociation 
contributions in the former samples. 



3 Data sets 

A comprehensive comparison has been carried out between recent HI and ZEUS measurements obtained 
with the three different methods. The data sets used are as follows. 1 

• Three data sets collected with the ZEUS detector in the years 1999 and 2000. Overlapping samples 
have been analysed with the ZEUS Leading Proton Spectrometer (termed "ZEUS LPS", based on 
a luminosity of 32.6 pb" 1 ) [16], with the LRG method ("ZEUS LRG", 62.4 pb~ x ) [16] and with 
the M x method, relying on the Forward Plug Calorimeter ("ZEUS FPC I", 4.2 pb" 1 [19] and 
"ZEUS FPC II", 63.4 pb^ 1 [20]). 

• A set of data collected with the HI Forward Proton Spectrometer ("HI FPS", 28.4 pb _1 ) [15] in 
the years 1999 and 2000. 

• A set of data collected with the HI detector in the years 1997, 1999 and 2000 and analysed with 
the LRG method ("HI LRG", 2.0 pb~\ 10.6 pb^ 1 and 61.6 pb" 1 for small, intermediate and 
large Q 2 , respectively) [9]. 

The HI LRG and FPS samples are statistically independent and are only weakly correlated through 
systematics. The three ZEUS samples also have different dominating systematics, but are not statistically 
independent. About 75% of events are common to both the ZEUS LRG and ZEUS FPC II data sets and 
35% of the ZEUS LPS events are also contained in the ZEUS LRG sample. 
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4 Proton dissociation background and corrections 

In proton dissociation processes at the lowest Mat, 
the dissociative system N often escapes entirely 
undetected into the forward beam-pipe. As Mat 
increases, it becomes more likely that dissocia- 
tion products are detected in the instrumentation 
most sensitive to forward energy flow. The LRG 
and Mx methods therefore do not distinguish low 
Mtv proton dissociation events from the case in 
which the proton is scattered elastically. Differ- 
ent cross-section definitions have been adopted, in 
which the proton dissociation contribution is ei- 
ther subtracted statistically, or else the quoted re- 
sults are integrated over a specific range of Mat. 
Since understanding the proton dissociation con- 
tributions and the corresponding corrections is 
fundamental to comparisons between the different 
measurements, a detailed discussion is presented 
in the following. 

In both the ZEUS LPS and the HI FPS 
analyses, the contribution from proton dissocia- 
tion events is negligible at small xjp <0.02. At 
the largest xn> values, it becomes kinematically 
possible for the detected leading proton to be the 

result of a decay of an N* or other proton excitation, the remaining decay products being unobserved. 
This background was estimated by ZEUS to contribute around 9% at xp = 0.1, using the Pythia Monte 
Carlo (MC) model [21]. In the HI FPS analysis, using the RAPGAP [22] implementation of the DlFFVM 
proton dissociation model [23], it was estimated to reach 3% at xjp = 0.08. 
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Fig. 2: (a) FPC energy and (b) LPS xl distributions for ZEUS 
proton dissociation samples (see text), with data compared to 
the tuned PYTHIA model, (c-e) Extracted fractions of proton 
dissociation events in the ZEUS LRG sample as a function of 
Q 2 , (3 and xtp after integration over the other variables [16]. 



'The comparisons here are restricted to published data and do not yet include the precise HI LRG and Mx method results 
obtained from 1999-2004 running [17]. 



Proton dissociation contributions in the LRG and Mx methods can be controlled using dedicated 
proton dissociation simulations tuned in M/v regions where dissociating protons leave signals in the de- 
tectors, and extrapolated into the Mjv regions where the dissociation products are typically not detected. 
In addition to this procedure, both HI and ZEUS use standard simulations of non-diffractive processes 
to control the small migrations of very high Mn or xp events into the measurement region, which occur 
due to inefficiencies of the forward detectors. 

4.1 ZEUS LRG 

In the recent ZEUS analysis, the Pythia simula- 
tion was tuned to proton dissociation signals. Two 
samples were selected by requiring activity either 
in the forward plug calorimeter (FPC) or at rel- 
atively low proton energy in the LPS. The sam- 
ples thus include the low M/v region in which 
proton dissociation products are invisible to the 
central detector. The generated distributions were 
reweighted in Mn, Mx and Q 2 to best describe 
the energy distribution in the FPC (£fpc)> an d the 
scattered proton energy fraction distribution (xl) 
in the LPS. Figures 2a and 2b show the compari- 
son of the reweighted Pythia model with the two 
proton dissociation samples as a function of these 
variables. Also shown in figures 2c-e is the result- 
ing estimate of the fraction of proton dissociation 
events in the LRG sample as a function of Q 2 , (5 
and xp. This fraction, obtained separately from 
the FPC and LPS samples, is constant at the level 
of 25%. 

The ratios of cross sections extracted from 
the ZEUS LPS and LRG data (the latter uncor- 
rected for proton dissociation background), are 
shown in figure 3. There is no significant dependence on Q 2 , ip or (3, illustrating the low xjp com- 
patibility between the two methods. The ratio averages to 0.76 ± 0.01(stat.)lg'g2(syst.)*g 05( norm O> 
the last error reflecting the normalisation uncertainty of the LPS data. The proton dissociation back- 
ground fraction in the LRG data is thus 24 ± l(stat.)^3(syst.)l 8 (norm.)%, in agreement with the result 
of the MC study, 25 ± l(stat.) ± 3(syst.)% (figure 2). Unless stated otherwise, the ZEUS LRG data 
are corrected by this factor in the following and thus correspond exclusively to the truly proton-elastic 
process. 

4.2 HI LRG 

The contribution from proton dissociation in the HI LRG analysis is constrained through the DlFFVM 
MC [23] model, normalised using the response to large Mjv events leaving signals in the forward and 
central detector components [9,24]. 

The data are corrected using DlFFVM to Mn < 1.6 GeV. The HI LRG data are then compared 
with the HI FPS measurement, in order to extract the proton dissociation cross section with Mn < 
1.6 GeV directly from the data. The ratio of the two measurements, after projection onto the Q 2 , ip 
and (3 axes, is shown in figure 4. There is no evidence for any dependence on any of the kinematic 
variables, as expected in the framework of proton vertex factorisation. The average value of the ratio is 
1.23 ± 0.03 (stat.) ± 0.16 (syst.), the largest uncertainty arising from the FPS efficiency. The result 
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Fig. 3: The ratio of the ZEUS LPS measurement (Mat = m p ) 
to the ZEUS LRG measurement before subtraction of proton 
dissociation background [16]. The lines represent the aver- 
age value of this ratio. An overall normalisation uncertainty 
of i g % is not included in the errors shown. 



is in good agreement with the DlFFVM estimate of 1.15_o' 8- The data and DlFFVM ratios translate 
into proton dissociation background fractions of 19 % and 13 %, respectively, consistent within the 
uncertainties. The similarity between the proton dissociation fractions in the raw HI and ZEUS LRG 
selections is to be expected given the similar forward detector acceptances of the two experiments. 
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Fig. 4: The ratio of the HI LRG measurement (corrected to Mn < 1.6 GeV) to the HI FPS measurement (Mn = m p ), after 
integration over the variables not shown in each case [15]. The lines represent a fit to the data assuming no dependence on any 
of the variables. An overall normalisation uncertainty of 13% is not included in the errors shown. 



4.3 ZEUS FPC 

The proton dissociation treatment is also critical in the Mx method, where the diffractive contribution 
is separated from the non-diffractive component in a fit to the inclusive \nM x distribution. Proton 
dissociation events with sufficiently large Mn for dissociation products to reach the FPC and central 
detectors lead to a reconstructed Mx value which is larger than the actual photon dissociation mass. 
The resulting distortion of the In M x distribution affects the diffractive contribution extracted in the fit 
if corrections are not made. According to the Sang MC model, the N system contaminates the Mx 
reconstruction for Mn > 2.3 GeV on average [25], and events in this Mn range are therefore subtracted 
using Sang before the \i\M x distribution is decomposed. The upper Mn cut in the SANG sample is 
denned by (Mn/W) 2 < 0.1, which leads to a variation of the subtracted fraction of events with W, the 
centre-of-mass energy of the photon-proton system. This contrasts with the LRG method, where MC 
studies confirm that the rapidity gap requirement efficiently eliminates proton dissociation at large Mn, 
the remaining fractional low Mjy contribution being independent of kinematics to good approximation 
(figures 3 and 4). 

Despite these difficulties, there is acceptable agreement between the ZEUS FPC data and the 
ZEUS LRG measurement. A global fit comparing the normalisations of the two data sets (after correcting 
the LRG data to Mn = m p ) yields a normalisation factor of 0.83 ± 0.04 to be applied to the ZEUS FPC 
results. This factor is comaptible with with expectations for the residual proton dissociation contribution 
based on the MC studies in sections 4.1 and 4.2. 

5 Cross section comparisons 

Due to their differing Mn coverages, the a,?® measurements from the different data sets are not directly 
comparable. However, assuming the factorisation of the Mn dependence which is suggested in the data, 
varying the Mn range should introduce only global normalisation differences, which can be estimated 
using the proton dissociation simulations. 
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Fig. 5: Comparison between the ZEUS Mx method ('FPC F and 'FPC IF) and ZEUS LRG method data [16]. As explained 
in the text, the Mx method data are scaled by a constant factor of 0.83 to account for proton dissociation contributions with 
M N < 2.3 GeV. 



5.1 Comparison between LRG and Mx methods 

ZEUS cross section measurements obtained with 
the LRG and Mx methods are compared in fig- 
ure 5. The LRG data are corrected to = M p °- 
as described in section 4. 1 and the relative nor- *~ ' 
malisation factor of 0.83 (section 4.3) is applied to 
the ZEUS FPC data to account for residual proton 
dissociation. The overall agreement between the , 
two measurements is good, apart from some dif- 
ferences at large x-p > 0.01. The Q 2 dependence 
of the Mx method data is also slightly weaker 
than that of the LRG data. 

5.2 Comparison between ZEUS LPS and HI 
FPS measurements 

The ZEUS LPS and HI FPS data are compared 
in figure 6. For this comparison, the ZEUS re- 



sults are extracted at the same (3 and Q values 
as HI and are therefore not affected by extrapola- 
tion uncertainties. The shape agreement is satis- 
factory and the overall normalisation discrepancy 
of around 10% lies within the large combined nor- 
malisation uncertainty of around 14%. 
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Fig. 6: Comparison between ZEUS LPS and HI FPS 
ments [16]. Normalisation uncertainties of ±10% 
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5.3 Comparison between ZEUS and HI LRG measurements 

The ZEUS LRG data are extracted at the HI and values, but at different Q 2 values. In order to 
match the M N < 1.6 GeV range of the HI data, a global factor of 0.91 ± 0.07, estimated with Pythia, 
is applied to the ZEUS LRG data in place of the correction to an elastic proton cross section. After this 
procedure, the ZEUS data remain higher than those of HI by 13% on average, as estimated with a global 
fit comparing the normalisations of the two data sets for Q 2 > 6 GeV 2 . This normalisation discrepancy 
is similar to that between the HI FPS and the ZEUS LPS data sets. It is in line with the errors due to 
the 8% uncertainty on the proton dissociation correction in the ZEUS LRG data and the 7% combined 
relative normalisation uncertainty between the two LRG data sets. 

In figure 7, the ZEUS results are scaled by a factor 0.91 x 0.87 (the factor 0.87 = 1-0.13 
normalising the ZEUS to the HI data) and compared with the HI LRG measurement. An excellent 
agreement between the Q 2 dependences is revealed throughout most of the phase space. There are small 
deviations between the dependences of the two measurements at the highest and lowest values. 
The results of the 'HI Fit B' NLO QCD DPDF fit to the HI LRG data [9] is also shown. It gives a good 
description of the data at large Q 2 . However, the extrapolation beyond the fitted region (Q 2 > 8.5 GeV 2 ) 
undershoots the precise new ZEUS low Q 2 LRG data, confirming the observation in [9] that a standard 
DGLAP fit to the lowest Q 2 data is problematic. 

6 A First Combination of Data Sets 

For easy future consumption at the LHC and elsewhere, it is desirable to combine the various HI and 
ZEUS diffractive DIS measurements into a single easily digestible HERA data set. Here we take the first 
steps towards this goal, by making a simple error-weighted average of the HI and ZEUS LRG data sets, 
ignoring correlations between the data points due to the systematic errors. LPS and Mx method data 
are not considered at this stage. For the purpose of this exercise, the ZEUS normalisation is fixed to that 
of HI as described in section 5.3 and shown in figure 7. The normalisation of the combined data thus 
has an uncertainty beyond the 10% level. Combinations can only meaningfully be made where there 
is basic agreement between the different measurements. Since this is not always the case at the lowest 
xtp values, we restrict the averaging to the xp = 0.003 and xp = 0.01 data. The combinations are 
performed throughout the measured Q 2 range, including the Q 2 < 8.5 GeV 2 region, beyond the range 
of the 'HI Fit B' parameterisation which is compared with the data. 

To account for the differences between the Q 2 binning choices, HI data points are adjusted to the 
ZEUS Q 2 values by applying small correction factors calculated using the 'HI Fit B' parameterisation. 
Where both collaborations then have measurements at a given (Q 2 , xjp,/3) point, a simple weighted av- 
erage is taken, using the quadratic sum of statistical and systematic uncertainties for each experiment, 
excluding normalisation uncertainties. At (Q 2 , xjp,/3) points where only one experiment has a measure- 
ment, that result is carried directly into the combined data set. 

The results of this averaging procedure 2 are shown in figure 8. They are indicative of the sort of 
precision which is achievable through combinations, with many data points having errors at the 3 — 4% 
level, excluding the normalisation uncertainty. At xjp = 0.01 the combined data agree well with the 'HI 
Fit B' DPDF parameterisation in its region of validity. At xjp = 0.003 the Q 2 dependences are also in 
good agreement with the parameterisation in the and Q 2 region of the fit, with the exception of the 
highest value, where the average is pulled towards the more precise ZEUS data. 

More sophisticated averaging methods may be used in the future, for example that [26] developed 
to perform similar combinations of inclusive HERA data, with a full systematic error treatment. No 
attempt has yet been made to extract DPDFs from the combined data. Based on the combined a®^ 1 and 
its Q 2 dependence shown here, no significant conflict is expected with the quark or gluon densities of 
existing DPDF sets such as 'HI Fit B' in the bulk of the phase space. However, small modifications are 

2 The authors take full responsibility for this combined data set; it is not an official Hl/ZEUS result. 
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Fig. 7: Comparison between the HI and ZEUS LRG measurements after correcting both data sets to Mm < 1.6 GeV and 
applying a further scale factor of 0.87 (corresponding to the average normalisation difference) to the ZEUS data. The measure- 
ments are compared with the results of the 'HI Fit B' DPDF extraction, which was based on the HI data shown. Further HI 
data at xp = 0.03 are not shown. 



likely to be necessary to the quark densities at small and large (3 values and understanding the region 
Q 2 < 8.5 GeV 2 remains a challenge. 

7 Summary 

HI and ZEUS diffractive DIS data obtained by various methods with very different systematics have 
been compared in detail. All measurements are broadly consistent in the shapes of the distributions. 
The comparisons between proton tagging and LRG method data internally to the two collaborations give 
compatible results on the proton dissociation contributions in the raw LRG selections. There is a global 
normalisation difference at the 13% level between the LRG measurements of the two experiments, which 
is a little beyond one standard deviation in the combined normalisation uncertainty. A similar difference 
is visible between the normalisations of the HI and ZEUS proton tagged data. 

A first step has been taken towards combining the two sets of LRG data, by arbitrarily fixing the 
normalisation to that of the HI data set and ignoring correlations within the systematic uncertainties in 
obtaining weighted averages. The results hint at the precision which might be obtained in the future with 
a more complete procedure. 
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Fig. 8: Combination of the HI and ZEUS LRG data following the procedure described in the text. The global normalisation is 
fixed to that of the HI measurement, in order most easily to compare the data with the 'HI Fit B' DPDF results. 
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